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We optimise the effective slcin-depth and resolution of Ag-Ti02, Ag-SrTiOs, and 
Ag-GaP multilayers for imaging with sub- wavelength resolution. In terms of trans- 
mission and resolution the optimised multilayers outperform simple designs based on 
combined use of effective medium theory, impedance matching and Fabry-Perot res- 
onances. For instance, an optimised Ag-GaP multilayer consisting of only 17 layers, 
operating at the wavelength of 490 nm and having a total thickness equal to one 
wavelength, combines 78% intensity transmission with a resolution of 60 nm. It is 
also shown that use of the effective medium theory leads to sub-optimal multilayer 
designs with respect to the trade-off between the skin depth and resolution already 
when the period of the structure is on the order of 40 nm or larger. 
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I. INTRODUCTION 



The presence of losses is a major limiting factor in the development of practical plas- 
monic devices. Here, we focus our interest on metal-dielectric multilayers for guiding optical 
wavefronts consisting of both propagating and evanescent planewaves which may carry infor- 
mation about objects with sub-wavelength dimensions^""— . The mechanism of such optical 
guidance may be attributed to coupling between surface plasmon polariton modes existing 
at the metal-dielectric interfaces, to negative refraction in photonic crystals, to effective op- 
tical anisotropy of the multilayer, to Fabry-Perot resonances, and typically to the interplay 
of these phenomena. In effect, an appropriately designed multilayer is capable of projecting 
the optical wavefront in between its external boundaries almost without diffraction includ- 
ing sub-wavelength details. Unfortunately, the permittivity of metals is complex-valued in 
the visible wavelength range resulting in a limited penetration depth inside the structure. 
For instance the skin depth of bulk silver is smaller than 20 nm in the range of visible 
wavelengths. 

There exist several ways of mitigating losses such as splitting the single-layer perfect 
lena^"— into a multilayer with thin metallic layers and compensating optical losses with 
gain,-'^^ using a high index dielectric with silver, which shifts the operational wavelength 
towards the red,-^^ or using semiconductors instead of silver for the wavelength ranges in 
between far-UV to far-IR.— Gain-tunable superresolution was recently considered along 
with the effect of gain on formation of optical vortices at the layer boundaries.- Reflection 
and transmission coefficients of the multilayer strongly depend on the termination condi- 
tions and in particular a symmetrical coating of the multilayer with a dielectric layer on 
both external boundaries leads to an increased transmission.- Simple yet successful designs 
of metal-dielectric superlenses for imaging with sub-diffraction resolution were designed with 
the effective medium theory (EMT)^'^ and recently refined beyond the second order Taylor 
expansion."^ The canalization regime of transmission- assumes that three effective permit- 
tivity and overall thickness conditions are fulfilled at the same time -the effective transverse 
permittivity of the multilayer is equal to that of the host medium Ex = 1, the axial permit- 
tivity is equal to = 0, and the total thickness of the multilayer satisfies the Fabry-Perot 
condition for transmission. However, the role of the impedance matching condition was later 
questioned.— In fact, multilayers with effective transverse permittivity 6^ = alone or the 
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FIG. 1. The geometry of a silver-dielectric multilayer. An elementary cell of the structure a consists 
of a silver layer with a thickness d2 and two dielectric coating layers, each with a thickness di/2. 

axial permittivity e^^ = alone were also proposed,- with larger absolute permittivity of 
metal than the absolute permittivity of dielectric in the first case and an opposite situation in 
the latter. Another approach assumes the use of complementary anisotropic slabs support- 
ing negative refraction, where each slab consists of a silver-dielectric multilayer.- Moreover, 
engineering of the point spread function of the multilayer with sub-wavelength full-width- 
at-half-maximum (FWHM) may be conveniently cast into the framework of Fourier optics 
adjusted to include the evanescent waves.-"— >^ Finally, numerical optimisation was used to 
improve FWHM in an impedance matched and impedance mismatched multilayer.— 

II. BACKGROUND 



The dispersion relation of a periodic multilayer for a TM-polarised monochromatic 
planewave takes the form^Si^ 

cos{kB a) = cos{kid\)cos{k2d2) — - { j-— + 7~~" ) sin{kidi)sin{k2d2), (1) 



2 \k2e1 kie2, 

where di and Si (with i = 1,2) denote the thickness and permittivity of the layers, a = di + d2 



is the period of the structure, ks is the Bloch wavevector, ki = \/k\ei — k^ is the component 
of the wavevector along the z-axis inside the i-th medium, and fco = 2vr/A is the free space 
wavenumber. A multilayer with N periods, total thickness L = N ■ a, and with symmetric 
termination is shown in Fig. [H An analogous dispersion relation for a TE-polarised wave was 
also presented by Wu et. al.,^^ however our interest is focused on metal-dielectric multilayers 
for the TM polarisation, since they enable SPP-enhanced transmission with sub-wavelength 
resolution. 
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When the layers are thin ki-di <^ 1, EMT makes it possible to approximate the structure 
with a uniaxially anisotropic slab with the effective permittivity tensor-i^'2^ 



ex 

, (2) 


where 

_ Eidi + £2C?2 _ di + d2 , s 

di + d2 £i di + 82 02 

are the effective permittivities for the directions parallel and normal to the layers surfaces, 
respectively. Equations and ([3]) may be derived from the electromagnetic boundary 
conditions by averaging electric fields E and D within a single elementary cell of the multi- 
layer.— 

In a lossless periodic layered structure may achieve an infinite value if the thicknesses 
di, d2, and permittivities £1, £2 satisfy the following relation:^ 

I -I- 

In a lossy structure Ez is always finite but may have a large magnitude which still enables to 
obtain a device that couples a broad spectrum of spatial frequencies, including both homoge- 
neous and evanescent waves, into a propagating mode. This mode is almost diffraction-free 
and in low-loss media it may be guided for a large distance. This enables the projection of 
sub- wavelength images from the front interface of the device onto the back interface and to 
obtain in-plane imaging with sub-wavelength resolution. 

Matching the impedances rj = \fjife of two media is a way to avoid reflections from their 
boundary. In a similar way, assuring in Eq. ([3]) that 1 or 

Exdx £2(^2 ^ c/i ^2 (5) 

makes it possible to eliminate reflections from the multilayer for normal incidence. A further 
increase of transmission and removal of reflections is possible for the Fabry-Perot condition 

fe^ ■ L m 



X 



-, where m = 1, 2 . . . (6) 



A 2 

Finally, when the size of a single dielectric layer is sufficient to form a cavity in between 
the metallic layers, cavity modes may be coupled enabling resonant transmission of a wave 



through the structure even when the metalhc layers are a lot thicker than the skin depth in 
bulk metal. The condition for a cavity mode to exist is 

A ^ 2 27r ' ^ ^ 

where I = 1,2 . . . and (fri2 is the phase of the complex reflection coefficient between materials 
with permittivities ei and £2 responsible for shortening of the cavity length. 

In order to measure the transmission efficiency of the multilayer we use the effective skin 
depth (5), defined as the distance at which the intensity of a normally incident wave decreases 
by a factor of the Euler constant e. Notably, sometimes the skin depth is defined with respect 
to amplitude rather than intensity, or the use of this term is restricted to the skin effect 
due to conductivity of metals, while here we use it in a broader sensed equivalently to the 
effective intensity decay rate and independently of its physical origin. Therefore the skin 
depth of a uniform or homogenised medium may be calculated using the simple formula, 

5 = A/47r3(n,), (8) 

where '^{ux) is the imaginary part of the effective refractive index in the direction parallel 
to the multilayer. For non-magnetic materials the effective index may be either defined 
using the effective permittivity tensor given in Eq. (E]) n^, = = ^/e^ or from the Bloch 

wavevector calculated using the dispersion relation = n^'-°^'^ = ksik^ = 0)/ko. We 

underline that the imaginary part of n^^°'^^ does not depend on the choice of the Brillouin 
zone and that the two definitions asymptomatically converge when a ^ A. 

As the measure of resolution we use the FWHM of the squared point spread function 
{\PSF\'^).— The point spread function is commonly used in Fourier optics to characterize 
the response of an optical system to a point source, whereas the \PSF\'^ is directly related 
to the intensity distribution in the image plane. 

III. OPTIMISATION RESULTS 

We focus on multilayers consisting of silver and three different dielectric or semiconduc- 
tor materials: Ti02, SrTiOs and GaP. Their respective dispersion characteristics used in 
modeling are taken from Johnson and Christy^E in case of silver and from the book by 
Palik2^ otherwise. The choice of these materials is based on their large permittivities in 
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FIG. 2. Skin depth of a) Ag-Ti02 and b) Ag~GaP periodic multilayer as a function of the fill 
factor of the structure. Particular curves correspond to different thickness of the elementary cell 
of the structure, varying from a = 40 nm to a = 200 nm. The black curve corresponds to the 
homogenised medium in accordance to the effective medium theory. The wavelength is equal to a) 
390 nm and b) 490 nm. The points marked at the curves refer to the values of the fill factor for 
which the cavity modes are supported in a single dielectric layer. Dashed lines indicate the filling 
fractions for diffraction- free propagation (in red) and for impedance matched to air (in black). 

the visible range and technical possibility of using them to prepare thin layers with several 
techniques.— i^""— 

In Fig.[2]we present the effective skin depth of Ag-Ti02 and Ag-GaP periodic multilayers 
as a function of the filling fraction of silver for several values of the period ranging from 
a = 40 nm to a = 200 nm and for the wavelengths of A = 390 nm and A = 490 nm, 
respectively. The skin depth for a periodic structure is calculated with Eq. ([8]) using the 
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effective index rix = obtained from the Bloch wavevector given from the dispersion 

relation (Eq. ([1])). Additionally, it is also calculated for an effective medium with rix = n^^'^, 
notably the two methods converge for a i— )■ 0. In the opposite limit, for large values of the 
period a, the multilayer behaves in a resonant way remaining transparent only for the size 
of dielectric layers corresponding to the cavity modes exited in the dielectric layers. The 
locations of these modes calculated with Eq. ([7]) are shown with circles. In a periodic 
structure with sufficiently thin metallic layers {dAg = a ■ (dAg/a) < 40 nm) these modes 
broaden into a band, which is reflected by the non-resonant behavior of skin depth in Fig. [2] 
for small fill-factors. Nonetheless, the specific shape of the field repeated in subsequent 
coupled cavities is also observed already within the band when the condition (JTj) applies. 
We will further demonstrate such a case for Ag — GaP and a = 60. Moreover, the filling 
factors which under EMT satisfy the conditions of diffraction- free propagation (Eq. (jlj)) 
and impedance matching (Eq. ([5])) with air are also marked in the same plot with vertical 
dashed lines. Further, we focus on structures with the period of a = 40 nm and a = 60 nm, 
primarily due to the technological feasibility of depositing thick (L ^ A) multilayers with 
such a period but also due to their increased skin depth, as compared with the multilayers 
with a smaller period (See Fig. [2]). 

Our goal in this paper is to optimize the resolution and skin depth of multilayers in the 
visible wavelength range. A periodic multilayer is defined by its total thickness L, period 
a, and the fill factor dAg/a- In the optimization procedure the resolution is measured for 
metal-dielectric stacks of the thickness equal to one wavelength L = X and the fill factor 
is optimized. The optimisation is of the trade-off character, since the effective skin depth 
decreases monotonically with the increase of the fill factor dAg/a while the FWHM of the 
\PSF\'^ takes minima only for specific values of the fill factor. We use the transfer matrix 
method (TMM) to calculate the PSF. In order to combine high resolution and transmission 
efficiency, we minimise the ratio FWHM/ 5 subject to the constraint that FWHM is less 
than A/2. 

In Fig.|3]we present the optimization results for Ag-Ti02, Ag-SrTiOs and Ag-GaP struc- 
tures, respectively. The resolution and the effective skin depth of the optimised multilayers 
are compared with those with the fill factor obtained using the EMT-based expression for 
diffraction- free propagation, Eq. (jl]). Clearly, for a broad range of wavelengths, optimization 
leads to the improvement of both resolution and transmission efficiency at the same time. 
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FIG. 3. Optimization results for Ag-Ti02 (a,b,c), Ag-GaP (d,e,f), and Ag-SrTiOa (g,h,i) multi- 
layers as a function of wavelength. Sub-figures include the fill factor (a,d,g), effective skin depth 
bjX (b,e,h), and a measure of resolution - FWHM of \PSF\'^ (c,f,i)- Fill factor is optimised for the 
sizes of period a = 40 nm and 60 nm (opt) and is compared to diffraction- free condition (EMT). 
Resolution is calculated for structures with the thickness of L = lA. 



8 



5 10 15 20 25 30 



FIG. 4. ID FDTD simulation with time-averaged Poynting vector inside the Ag-Ti02 structure for 
o = 40 nm (lower blue line) and for a = 60 nm (upper red line) at the wavelength of A = 390 nm. 
The fill factors of both structures are optimised. The intensity level of 1/e is marked with a dashed 
line. The inset shows a magnified part of the plot. 

The largest skin depth is achieved for Ag-Ti02 structure for A = 390 nm, and for Ag-GaP 
structure for A = 490 nm. They are equal to 5 = 14 A and to 5 = 6.5 A, respectively. In 
both cases the resolution is on the order of 0.12 — 0.15 A and it is little affected by a change 
in a. For Ag-SrTiOa the skin depth does not exceed 1.3 A for any wavelength, which makes 
strontium titanate less interesting for fabricating thick multilayers, although the maximal 
value of skin depth corresponds to the resolution better than 0.1 A. The skin depth is larger 
for structures with the larger period a = 60 nm than for those with a = 40 nm. We illus- 
trate this increase with a one-dimensional FDTD simulation for the Ag-Ti02 stack shown 
in Fig. m The plot includes the time-averaged Poynting vector inside an infinite structure. 
For both periods, light decays exponentially, however the larger thickness of elementary cell 
results in a slower decay rate. The inset in Fig. H] contains a magnified part of the same 
plot showing that the energy is absorbed only in the silver layers, since titanium dioxide is 
practically lossless in the visible range. The decay rate obtained with FDTD is in perfect 
agreement with the value of effective skin depths calculated with the dispersion relation ([T]). 

Good insight into the transmission mechanism involved in the imaging through the mul- 
tilayer is gained from the internal field distribution. With this aim, in Fig. [5] we summarise 
the field profiles for several designs of Ag-GaP structure at A = 490 nm with the total 
thickness of L = lA, either with a = 60 nm or with extremely thin layers a 0. We 
compare the following structures: the one which we have optimised, the one consisting of 




FIG. 5. Transmission of a normally incident plane wave through Ag-GaP multilayers with the 
total thickness L = A at the wavelength of A = 490 nm. a) multilayer optimised for a = 60 nm; 
b) coupled-cavity multilayer with a = 60 nm; c) homogenised version of (a); d) multilayer with 
maximised £z] e-h) multilayers satisfying the Fabry-Perot condition of the order m = 1,2,3,4. 
Multilayer with m = 2 is at the same time impedance matched to air. 
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TABLE I. Comparison of the effective skin-depths {6), resolutions {FWHM{\PSF\'^)) and inten- 
sity transmission coefficients for normal incidence (T) of the Ag-GaP multilayers with the total 
thickness L = A at the wavelength of A = 490 nm. 



Description of the multilayer 


rl . In 

aAg/a 


^ 1 \ FWHM 
^1^ A 


T 
1 


Optimised for a = 60 nm 


0.270 


6.31 0.12 


0.78 


Coupled cavities with a = 60 nm 


0.522 


2.00 0.58 


0.61 


Optimised and homogenised 


0.270 


4.22 0.76 


0.42 


Diffraction- free 


0.412 


2.23 0.08 


0.58 


Fabry-Perot m = 1 


0.576 


0.45 0.08 


0.10 


Impedance-matched, F-P m = 2 


0.542 


0.90 0.11 


0.33 


Fabry-Perot m = 3 


0.487 


1.47 0.11 


0.49 


Fabry-Perot m = 4 


0.408 


2.27 0.08 


0.59 



coupled-cavities (with the width of ^2 defined using Eq. for the fundamental mode), the 
homogenised analogue of the optimised one (split into more layers with the filling fraction 
and total thickness fixed), a diffraction- free structure (fill factor calculated with Eq. @), 
and Fabry-Perot etalons of the order m = 1, 2, 3, 4 (with m = 2 also satisfying the condition 
for impedance matching ([5])). The corresponding filling fractions, skin depths, resolutions 
and intensity transmission coefficients are compared in Tab. HI The immediate conclusion 
from this summary is that the optimised multilayer does not resemble the theoretical designs 
neither in terms of the filling fraction nor in terms of the internal field distribution. This 
is a strong argument for using numerical optimisation in place of simplified models based 
on EMT when it is important to combine a good transmission efficiency, a subwavelength 
resolution and to keep the total number of layers technologically feasible. 

In Fig. |6] we present the amplitude and phase of amplitude transfer functions (TF) for 
Ag-Ti02 multilayer (at A = 390 nm) and for Ag-GaP multilayer (at A = 490 nm). The 
period a and the total thickness of structure L are fixed at 60 nm and 3 A, respectively. 
Horizontal cross-sections of Fig. [H] include the TF calculated for a range of fill factors, and 
these TF depend on k^/ko. 

For propagating waves k^/ k^ is the sine of the angle of incidence. The values of wavevector 
kx/ko > 1 refer to the evanescent waves in air. Optimal fill factors for both structures are 
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FIG. 6. Dependence of the amplitude transfer function TF (in horizontal cross-sections of the 
figure) on the fill factor: a, c) logarithm of amplitude, and b, d) phase. Two structures were 
considered: a, b) Ag-Ti02 multilayer (for A = 390 nm) and c, d) Ag-GaP multilayer (for A = 
490 nm). The total thickness of both structures is equal to 3 A and the thickness of the elementary 
cell a = 60 nm. The optimal fill factors are 0.22 for Ag-Ti02 and 0.28 for Ag-GaP. 

0.22 and 0.28 for Ag-Ti02 and Ag-GaP multilayers, respectively. At the same time, for 
these values of the fill factor, the phase of TF in Fig. [6] is the flattest, while the amplitude is 
slowly varying. The phase is practically independent of kx/ko and is similar for propagating 
and for evanescent components of the spatial spectrum enabling diffraction-free propagation 
with sub-wavelength resolution. Moreover, sub-wavelength resolution may be obtained even 
if the overall thickness of the multilayer is much larger than the wavelength. 

Let us now illustrate the performance of the optimised multilayer. Up to now we have 
neglected any possible interaction between the source of the sub-wavelength-sized wavefront 
and the multilayer. In practice, the source may be realised using a metallic mask with sub- 
wavelength aperture, a kinoform with sub-wavelength relief, a plasmonic waveguide coupled 
to the multilayer, or with a SNOM probe. In any case, the source needs to be put in the 
near field or preferably be attached to the multilayer, therefore causing multiple reflections 
in between both elements. The system consisting of the source coupled to the multilayer 
analyzed as a whole, may be no longer shift invariant, although idealized source models such 
as a soft (uncoupled) source, or hard (strongly coupled) magnetic and electric sources may 
still be included in the same framework of LSI — . The three source models become equivalent 
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FIG. 7. FDTD simulation of Ag-GaP multilayer consisting of 25 elementary cells with a = 60 nm 
and fill factor 0.28. The structure is illuminated with a plane wave with A = 490 nm diffracted at 
a subwavelength aperture with the diameter equal to A/20, a) The z-component Sz of Poynting 
vector; b) intensity profiles in three cross-sections: at input and output interfaces of the multilayer 
and at z = IX from the input interface. 

in the absence of reflectiona^i. Even though, a soft source model is assumed by us in this 
paper, in case of the optimised multilayer, thanks to small reflections, the FWHM of PSF 
is nearly the same for the three source models. 

In Fig. [7] we present the results of a two-dimensional FDTD simulation of a Ag-GaP 
structure consisting of 25 elementary cells with a = 60 nm and fill factor equal to d-Ag/cL = 
0.28. We show the longitudinal component of the Poynting vector and intensity profiles at 
three cross-sections perpendicular to the propagation direction. The structure is illuminated 
with a monochromatic planewave with A = 490 nm diffracted at a sub-wavelength aperture 
with the diameter equal to A/20. The aperture is made in a perfect metal screen having the 
thickness of 2 nm located at the distance of 10 nm from the multilayer. Transmission through 
the aperture is not enhanced by resonant effects. Resonant transmission would certainly 
complicate our analysis, whilst our main focus here is on the transmission properties of the 
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multilayer rather than on the efficient wavefront modulation at the scale of 20 nm. This way, 
the presented results are general in the sense that they are little affected by the choice of 
metal used to make the mask, provided that its thickness is sufficient to make it opaque. For 
instance, a chromium mask with the thickness of 50 nm could be used instead the perfect 
conductor giving similar results. 

The sub-wavelength dimension and shape of the wavefront is preserved within the struc- 
ture despite the continuous intensity decay due to the losses in metallic layers. Moreover, 
the FWHM at the distance of lA is consistent with the optimisation results obtained with 
the transfer matrix method and presented in Fig. [Sf. 

IV. CONCLUSIONS 

We have optimised the effective skin-depth and resolution of Ag-Ti02, Ag-SrTiOa, and 
Ag-GaP multilayers for imaging with sub-wavelength resolution. We have shown that mul- 
tilayers with a period of a = 40 nm 8% ■ A) or larger, designed using the effective medium 
theory (Eqs. dl]) , ([5]) , ([6]) ) are suboptimal in terms of both resolution and effective skin depth 
at the same time. Still, multilayers with a large thickness of layers and with silver lay- 
ers slightly thicker than the skin depth may combine a good trade-off between resolution 
and transmission. Optimal multilayer designs with a period of a = 40 nm or a = 60 nm 
outperform multilayers based on combined use of effective medium theory (with a/A t— t- 0, 
N = L/a and L = const), impedance matching and Fabry-Perot resonances in terms of 
transmission, with a certain deterioration of resolution. For instance, an optimised Ag-GaP 
multilayer consisting of only 17 layers, operating at the wavelength of 490 nm and having 
the total thickness equal to the wavelength, combines the intensity transmission coefficient 
of 78% with the resolution of 60 nm. Using EMT, the resolution may be improved to 37 nm 
(See Tab. [T]) but to achieve this it is necessary to increase the number of layers several- fold 
which is technologically challenging. We have shown that an optimised structure allows for 
diffraction-free guidance of a sub-wavelength sized beam of light for at least L = 3 ■ A and 
that the transmission is not directly based on Fabry-Perot resonances. At the same time, 
the optimal design does not resemble the EMT-based multilayer nor the coupled-cavity mul- 
tilayer in terms of either the ffiling fraction or the internal distribution of field making an 
argument for the use of numerical optimisation instead of relaying on simplified theoretical 
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designs. 
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